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Abstract 

We present a consistent theory of superconductive tunneling in single-mode 
junctions within a scattering formulation of Bogoliubov-de Gennes quantum 
mechanics. Both dc Josephson effect and dc quasiparticle transport in volt- 
age biased junctions are considered. Elastic quasiparticle scattering by the 
junction determines equilibrium Josephson current. We discuss the origin 
of Andreev bound states in tunnel junctions and their role in equilibrium 
Josephson transport. In contrast, quasiparticle tunneling in voltage biased 
junctions is determined by inelastic scattering. We derive a general expres- 
sion for inelastic scattering amplitudes and calculate the quasiparticle current 
at all voltages with emphasis on a discussion of the properties of subgap tunnel 

current and the nature of subharmonic gap structure. 
PACS: 74.50.+r, 74.20.Fg, 74.80.Fp 
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I. INTRODUCTION. 



The tunnel Hamiltonian model |l| has for many years been a main theoretical tool 
for investigation of tunneling phenomena in superconductors 0. However, interpretation 
of recent experiments on transmissive tunnel junctions and complex superconductor- 

semiconductor structures |],[7| require more detailed knowledge of the mechanisms of the 
superconductive tunneling than the tunnel model is able to provide. Particularly informative 
are experiments on superconducting quantum point contacts with controlled number of 
transport modes and transparency, such as controllable superconducting break junctions || 
and gate controlled superconductor-semiconductor devices ||. Since only a few transport 
modes with controlled transparency are involved in the tunnel transport, the experiments 
provide precise and detailed information which can be directly compared with theory. 

The first attempts to develop a theory of superconductive tunneling beyond the tunnel 
Hamiltonian model fl0|-|13lj were made in generalization of methods applied to SNS junc- 
tions [M| ) 15 and superconducting constrictions fl6| , [i7H based on Green's function techniques. 
In these theories, the junction Green's functions are directly found from Green's function 
equations supplemented by special boundary conditions representing the tunnel barrier, or 
by matching the superconductor and insulator Green's functions at the superconductor- 
insulator boundaries. 



In the first works on the Josephson effect in SNS junctions Hl8l , |T9f another way of cal- 
culation has been used, based on expansion over eigenstates of the Bogoliubov-de Gennes 
(BdG) equation p0[ . A similar method has been also applied to SIS tunnel junctions |21 



and superconductor-semiconductor junctions ||22|| . In the absence of inelastic scattering the 
method of using the BdG equation gives the same results as the Green's function technique 
pTfl . One might then expect that the Josephson effects in superconducting junctions can be 
explained on a rather simple quantum mechanical level. Following this idea, the quantum 
mechanical approach has been successfully applied to calculation of the dc Josephson current 
in different kinds of mesoscopic weak links [fM|-|2"8|| and tunnel junctions |29|-|31|1 . The first 
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application of the method to voltage biased junctions was done by Blonder, Tinkham, and 
Klapwijk who considered quasiparticle tunneling in SIN junctions as a scattering problem in 
BdG quantum mechanics . Later on, the quantum mechanical approach has been found 
to be helpful in investigations of more complex phenomena of quasiparticle transport and ac 
Josephson effect in voltage biased SNS junctions P2| , mesoscopic SIS tunnel junctions |5B 
and mesoscopic constrictions |34| . 



The quantum mechanical approach based on the BdG equation is adequate for describing 
the physical situation in mesoscopic junctions, where inelastic scattering effects are weak 
and coherent electron dynamics is of main importance. Moreover, due to the effect of 
quantization of transverse electron modes in mesoscopic junctions f24j , |35|l , ID models for the 
current transport through the junction may be appropriate. 

In this paper we present a consistent quantum mechanical theory of superconductive 
tunneling in a one-mode quantum constriction, Fig.l. We consider the dc Josephson effect 
and also dc quasiparticle tunneling in voltage biased junctions. In the latter case we focus 
attention on a detailed calculation of the subharmonic gap structure (SGS) of the tunnel 



current IBS 



Following the Landauer approach pEfl , we consider superconducting electrodes as equi- 
librium reservoirs which emit quasiparticles into the constriction. Scattering by the junction 
goes into two channels: (i) the normal channel with outgoing quasiparticles remaining in 
the same branch of the quasiparticle spectrum, and (ii) the Andreev channel where the 
quasiparticles change branch due to electron-hole conversions. The current in such a picture 
results from the imbalance of currents carried by scattering states originating from the left 
and the right reservoirs, the magnitude of the current being proportional to the transmission 
coefficient D of the tunnel barrier. 

The imbalance of currents in superconducting junctions can be created in two ways: by 
establishing a difference of the phases of the order parameters in the left and right electrodes, 
or by applying a voltage bias. The basic fact concerning the flow of equilibrium current in 
the presence of a phase difference, established by Furusaki and Tsukada p{| , is that a bulk 



supercurrent, when approaching the tunnel interface, is transformed into current flowing 
through superconducting bound states which appear at the tunnel interface in the presence 
of this phase difference [[J7[] and which provide transmission of the Cooper pairs through the 
tunnel barrier. The balance among currents of different scattering states is not violated, 
although the scattering amplitudes depend strongly on the phase difference. 

Application of a voltage bias gives rise to more far reaching consequences than just im- 
balance of the elastic scattering modes: the scattering states themselves are modified in 
a non-trivial way. This follows from the fact that the scattering amplitudes, being phase- 
dependent in equilibrium, become time-dependent in accordance with the Josephson relation 
|38|, d(p/dt = 2eV, when voltage is applied. Thus, in the presence of a dc voltage the su- 
perconducting junction behaves as an effective nonstationary scatterer whose transmissivity 
oscillates. This property of superconducting junctions gives rise to ac Josephson effect; how- 
ever, it is also significant for dc quasiparticle transport because the quasiparticle transmission 
through such a scatterer is inelastic. 

The physical mechanism of inelastic quasiparticle transmission through voltage biased 
superconducting junctions has been first considered in SNS junctions where it has been 



explained in terms of multiple Andreev reflections (MAR): the normal quasiparticles con- 
fined between superconducting walls are permanently accelerated by the static electric field 
due to sequential electron-hole conversions at the NS interfaces, similarly to acceleration 
of the electrons in an ordinary potential well by a time dependent electric field. Similar 
arguments can be extended to the tunnel junctions |39|; however, in tunnel junctions the 



scattering theory approach is more appropriate because of the quantum nature of quasi- 
particle transmission through the atomic-size tunnel barrier. This introduces a side-band 
spectrum of scattered waves where side-band energies are shifted with respect to the energy 
of the incident wave by integer number of quanta of the scatterer frequency ||33|| . Such an 
approach is familiar in the theory of quantum scattering by oscillating potential barriers in 
normal tunnel junctions (see e.g., Refs. pp|,^T| and references therein). 



The tunneling through all the inelastic channels (normal and Andreev) constitutes a 



complete picture of superconductive tunneling in biased Josephson junctions - the incoherent 
part of side band currents corresponding to the dc quasiparticle current and the side band 
interference currents corresponding to the ac Josephson current. An important aspect of 
this picture is that the Andreev bound states are involved in the current transport together 
with the extended side band states, giving a multiparticle character to the superconductive 
tunneling in the subgap voltage region. This multiparticle origin of the subgap tunnel current 
was first pointed out by Schrieffer and Wilkins f42j. 

The structure of the paper is the following. After formulation of the problem and dis- 
cussion of the quasiclassical approximation in Sec. II, we consider the problem of elastic 
scattering in Sec. Ill as a starting point for construction of inelastic scattering states in 
biased junctions. The solution of the elastic scattering problem allows us to calculate dc 
Josephson current, which is done for completeness in Sec. IV. In Sec. V we construct 
inelastic scattering states and derive a continued-fraction representation for the scattering 
amplitudes. Sec. VI is devoted to derivation of the nonequilibrium current. In Sec. VII 
we discuss the origin of the excess tunnel current in the large bias limit. In Sec. VIII we 
present a general analysis of the subgap tunnel current. Finally, the SGS is analyzed in Sec. 
IX. 



II. FORMULATION OF MODEL. 



We consider a superconducting quantum constriction with adiabatic geometry f43fl : the 
cross section varies smoothly with the coordinate x on the scale of the Fermi electron wave 
length, 1/pf, and the size of cross section is comparable with the Fermi electron wave length, 
Fig. 1. The length L of the constriction is assumed to be smaller than the superconducting 
coherence length £ , 

I/PF < L < fo- (2.1) 
The Hamiltonian of the constriction is assumed to have the form: 
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H 



(p-a z eA(r,t)y 



2m 



+ U(r)- n 



o z + [V(x) + e<p(r, t)} a z + A(f, t), 



(2.2) 



where U(r) is the potential confining electrons within the constriction, V(x) is the potential 
of the tunnel barrier, A{r, t) and <f{r,t) are electromagnetic potentials, and A(f,t) is the 
off-diagonal superconducting order parameter given by the matrix: 



A 



Ae^/ 2 
Ae- ,; */ 2 



(2.3) 



We assume the junction to be symmetric. The choice of the units corresponds to c = % — 1. 

It is convenient to eliminate the phase of the superconducting order parameter x(r*, t) in 
Eq. (|2.3|) by means of a gauge transformation: 



(2.4) 



which allows us to introduce a gauge invariant superfluid momentum, p s = Vx/2 — eA and 
an electric potential $ = x/2 + e<y9. 

There are different scales of change of potentials in Eq. ( |2.2j ): one is an atomic scale over 
which the confining potential U(r±) and the potential of the tunnel barrier V(x) change. 
Other scales are related to change of superconducting order parameter, electromagnetic field 
penetration lengths and length of the contact: all these lengths are large in comparison with 
the atomic length. It is convenient to separate these two scales, introducing quasiclassical 
wave functions JPJ which vary slowly on an atomic scale, and including rapidly varying 
potentials in a boundary condition for quasiclassical wave functions. To this end we assume 
the solution \l/(r*, t) of the Bogoliubov-de Gennes equation |20 



with the Hamiltonian of Eq. Q2.2|) , to have the quasiclassical form 



(2.5) 



(2.6) 
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where ip± is the normalized wave function of the quantized transverse electron motion with 
the energy E±, 

— * 

(— + V(r±,x))tjj± = E±(x)ijj± } ip±(r± = oo,x) = 0, 
2m 



and p is the longitudinal momentum of the quasiclassical electron, p(x) = J 2m(// — E±(x)); 
(3 = ± indicates the direction of the electron motion. We assume here that the constriction 
has only one transport mode; an extension to the case of several unmixed modes consists 
of additional summation over all transport modes in the equation for the current. The 
coefficients ip? in Eq. ( |2.6| ) describe the wave functions slowly varying in the x direction and 
satisfy the reduced BdG equation 



(2.7) 



in the the left (L) and the right (R) electrodes; v = p/m. The potentials p s and $ describe 
the distributions of electromagnetic field and supercurrent in the electrodes. In the point 
contact geometry these quantities are small due to the effect of spreading out of current 
16|)45l) and we will omit them, p s = $ = 0. For the same reason, deviation of spatial 



distribution of the module of the order parameter A from constant magnitude is small in 
the point contacts, and we will neglect it, A = const. 

The functions ip^ R are matched at the constriction by the boundary condition |HJ (see 
also Appendix A): 



(2.8) 





= V 








cltj X 











with a matching matrix V 



V 



r de ia ^ 12 



(2.9) 



d and r are the normal electron transmission and reflection amplitudes due to the barrier, and 
is a gauge invariant difference of the superconducting phases of right and left electrodes: 
<p = X-r(O) — The matching matrix in Eq. (|2.9| ) satisfies the unitarity condition 



VV f = I. 



(2.10) 



The boundary condition in Eqs. ( p78|) , fl2.9|) is analogous to the boundary condition used 
in the quasiclassical Green's functions techniques (see e.g., pd] , p6| ). This is the simplest 
equation for coupling of superconducting electrodes, while retaining the main features of 



the Josephson effect, except for effects of the resonant tunneling [30,47,48 



III. ELASTIC SCATTERING 

In the absence of time dependence in the phase difference at the junction, = 0, Eqs. 
( f2~7D and ( |2.8p describe elastic scattering of quasiparticles. The scattering states are to be 
constructed using stationary solutions of Eq.( f2~7l) which correspond to elementary propa- 
gating waves with energy \E\ > A: 

/ e <57/2 \ 



u 



J 



y/2 cosh 7 
where 

\E\ -\- £ 

£ = VE 2 - A 2 , e< = 1 \ — , a = signE, a = ±, 5 = aa 



(3.1b) 



(3.2) 



The vector function ue is normalized, (u , u) = 1, the brackets meaning scalar product in 
electron-hole space. In Eq. ( |3.1| ) there are four elementary waves which correspond to the 
same energy, as illustrated in Fig. 2, which are labeled by quantum numbers (3 (direction of 
the Fermi electron momentum) and a = sign(\p\ — pp) (electron or hole-like branch of the 
quasiparticle spectrum). The direction of propagation of each elementary wave is determined 
by the sign of the probability current. The probability current density j p , defined by the 
conservation law (continuity equation) d\ip\ 2 /dt + dj p /dx = for the BdG equation Eq. 
( pT7|) , has the form j p = (ip,a z ip). For the elementary waves in Eq. (|3.1|) one obtains the 



S 



explicit result j p = /35tanb.7. According to this formula, the relation S = f3 is met for for 
the waves propagating from left to right, 5 = —f3 for the waves propagating from right to 
left. Therefore the incoming waves from the left (L) and right (R) have the form 



L : e ia ^ /v)x u E) R : e~ i(T ^ lv)x u^ \ 
while the outgoing waves have the form 

L : e- ia ^ /v)x UE P , R : e ia ^ /v)x u E . 



(3.3) 



(3.4) 



Correspondingly, the incoming quasiparticle can be scattered into four outgoing states: two 
forward scattering states and two back scattering states. One of the reflected waves belongs 
to the same (electron- or hole-like) branch of the quasiparticle spectrum as the incoming 
wave and constitutes the normal scattering channel, while the other reflected wave changes 
spectrum branch and constitutes the Andreev channel. In a similar way transmitted waves 
constitute normal and Andreev channels. The structure of the scattering states then becomes 



V ^ 2 J 



e i*(Z/v)x u - + 



U 



(3.5a) 







<8 \ 









e a E 



(3.5b) 



(for brevity we have omitted the time dependent factors e~ lEt ). In Eqs. ( |3.5| ) the index 
j = 1(2) corresponds to a hole-like quasiparticle coming from the left (right), while index 
j = 3(4) corresponds to an electron-like quasiparticle coming from the left (right). According 
to the structure of the matching matrix Eq. ( |2.9| ), there is the following symmetry between 
the scattering states j = 1 and j = 2: 



(3.6) 



















w = 




(-0), 




(</>)= 




to 


K a J 


i 
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The analogous symmetry exists also for the scattering states j = 3, 4. Using the unitarity of 
the matching matrix Eq. ( |2.9| ) one can find the following relation between scattering states 
j = 3 and j — 1 : 



(7, r, d) 



(A 



(-l.r'.d"), 



(A 



(l, r, d) 



V 6 A 



(-7,r*,d*). (3.7) 



These symmetry relations allow us to find all the scattering amplitudes if one of the scattering 
states is known. 

Let us find the explicit scattering amplitudes for the scattering state j = 1. After 



substituting Eqs. ( |3.5|) into Eq. (|278|) it is convenient to split the resulting equation, using 
the orthogonality condition, (u + , o z u~) = 0, into two independent equations for the normal 
scattering amplitudes c, / and for the Andreev scattering amplitudes a, b, 



(A 



\ / 



(u ,cr z Vu ) 



(3.8a) 



(u + ,a z u + ) 



(u + , a 7 Vu 



\0> 



(3.8b) 



Calculating the scalar products in Eqs. ( |3.8| ) we find explicit expression for the Andreev 
amplitudes in terms of the normal ones, 



^ a\ idsm((f)/2) 



sinh7 



The solution of the first equation in Eq. ( |3.8|) is given by 



Cl 



r sinh 2 7 



Ji 



fi 



c?sinh7sinh(7 + i<p/2) 



(3.9) 



(3.10) 



where 



d 



Z = (R sinh 2 7 + D sinh(7 + i<j)/2) sinh(7 - z0/2) 

d* v 



(3.11) 
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D = \d\ 2 is the normal electron transmission coefficient of the tunnel junction, and R = 
\r\ 2 = 1 — D is the normal electron reflection coefficient. It follows from Eqs. (|3.9| ), ( |3.10| ) 
that if there is no phase difference across the junction, = 0, the Andreev scattering channel 
is closed: a = b = 0. It is worth mentioning that the Andreev reflection is also absent if 
the normal transparency of the junction is equal to zero, D = 0. If, on the other hand, the 
junction is completely transparent for normal electrons, D = 1, there is no Andreev forward 
scattering, b = c = 0. 

In the presence of a phase difference at the junction the quasiparticle scattering is ac- 
companied by the appearance of superconducting bound states ||37|| . One can establish the 
existence of bound states by investigating the poles of the scattering amplitudes, Eq. ( [3.10D , 



at imaginary 7 corresponding to energies lying inside the gap \E\ < A. Assuming 7 — > 27 
in Eq. ( |3.11|) , we have the dispersion equation Z{i^) = 0, or 



sin 2 7 = D sin 2 0/2. (3.12) 

The bound states correspond to a positive magnitude of sin 7: A sin 7 = Im £ > 0. This 
condition selects two roots: 



7 = 70 = arccos(v D sin 0/2), 7 = 7r — 70, (3.13) 

or 

E(4>) = ±Ay/l -£>sin 2 0/2. (3.14) 

The wave functions of the bound states are to be constructed from elementary solutions of 
Eq. (|2.7|) with \E\ < A which decay at x = ±00, 

= e- iEt - <x/v W E , (3.15a) 

= e-^ +( >/^, (3.15b) 

where 
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( p iwy/2 \ 



u; 



V2 



ae 



—ivf/2 



\E\+K 



C = VA 2 - E 2 , v = (3a. 



(3.16) 



The bound state ansatz has a form similar to the outgoing part of the scattering states 
Eq. Q3.5|) with the coefficients satisfying the homogeneous equations of Eq. ( |3.8|) . These 
coefficients are: 



(isin(7 + 0/2) 
/ = : C, 



r sin 7 



/ \ ( 
1 a 1 dsin(0/2) 



sm7 



./ 

\- C J 



(3.17a) 



(3.17b) 



with 7 given by Eq. Q3.12 ). We note that the bound state spectrum is nondegenerate. The 
coefficient c in Eqs. ( 3.17 ) is obtained from the normalization condition for the bound state 
wave function, 

/ dV r dx\V\ 2 = -A\a\ 2 + \b\ 2 + |c| 2 + |/| 2 ) = 1, 

J J — oo (, 

which yields 

J Dsin 2 (7 + 0/2)\ _1 



|c| = A sin 7 1 + 



R 



sin 2 7 



(3.18) 



What is the origin of the bound states in a tunnel junction? According to Eq. 
one can regard these states as resulting from hybridization of the bound states in the short 
ballistic constriction [E^J] due to the normal electron reflection by the barrier (cf. effect of 
impurities in the SNS junction J2RE_B|). Let us consider a smooth constriction with the length 
exceeding the coherence length, L>(o. In such a constriction the supercurrent density and 
the superffuid momentum are related by the local equation, J s (x) = (e/m)N s p s (x), and 
they are both enhanced in the neck of the constriction due to current concentration (for 
simplicity we neglect the effect of suppression of the superfluid electron density N s by the 
supercurrent). The local quasiparticle spectrum in the presence of supercurrent has an 
additional contribution ±vfP s (x) p0| , which gives rise to a shift of the local energy gap, 



12 



Fig. 3. The spatial bending of the gap edges forms the potential wells at E < 0(E > 
0) for quasiparticles with electron velocities directed along (opposite) the current. The 
bound states in these potential wells are similar to the Andreev bound states in the SNS 



junctions j29|, the difference being that here the bound states are caused by the spatial 
nonhomogeneity of the phase of the order parameter, while the original Andreev states 
are caused by the spatial nonhomogeneity of the modulus of the order parameter. With 
decreasing length of the constriction the number of the bound states in the well decreases. 
The short Josephson constriction corresponds to an infinitely narrow and deep 5-potential 
well which contains only one Andreev level ||24|| . 



IV. DC JOSEPHSON CURRENT 



A convenient expression for the tunnel current results from statistical averaging of the 



current operator written in the Nambu representation p0|: 



I(x,t) = — \ (P-P') I d 2 r± \6{r-f*) -Tr (tf(f,t) ^(f^t)) 



(4.1) 



where ^ is a two-component field operator: 

/ 



V 



/ 



(4.2) 



and Tr is a trace in electron-hole space. The angular brackets in Eq. ( ^4.1| ) denote a thermal 
average of the one-particle density matrix of the superconductor ||51|| . In equilibrium this 
matrix has the form 



(*(f) *t(f*)) =£^(r>H-£A)*l' 



(4.3) 



where ^f\(r) are the eigenstates of the stationary BdG equation Eq. ( |2.5| ) with the quantum 
numbers A. We note that the definition of Fermi distribution function rtp here corresponds 
to the distribution of holes in the normal metal: in the ground state all energy levels above 
the Fermi level (E > 0) are occupied, while energy levels below the Fermi level (E < 0) are 



13 



empty (see also the discussion in the next section). In the quasiclassical approximation, Eq. 
the average tunnel current calculated at the middle of the junction has the form 



a a 



(4.4) 



The current in Eq. Q4.4| ) can be calculated either at the left or the right side of the junction 



because, due to the unitarity of the matching matrix V in Eq. ( p,10|) , the equality 



k#l 2 -|^l 2 H^I 2 -|V> 



R\ 



(4.5) 

holds for the each eigenstate. The current in Eq. ( |4.4|) consists of contributions from both 
the scattering states and the bound states 



dE\E\ 
e\>a 2n£ 



n F {-E)Y,!A E ) ~ E n F (-E)I bound (E), 



(4.6) 



\E\<A 



i{E) = eY J m p m 2 - 

When calculating the contribution from the scattering states, it is convenient to consider 



the transmitted current of each scattering mode: 



e(N 2 "l//) J = 1,3 
e(M 2 -hf) J =2,4 

The symmetry relations Eqs. (|3.6|), ( |3.7| ) yield 



(4.7) 



h(E) = h(E), I 2 (E) = I 3 (E), I 2 {E) = -h{E). 



(4.8) 



Thus the currents of all the scattering states with a given energy cancel each other in 
equilibrium fl52|. Substitution of equations Eq. ( |3.17|) , (|3.18|) into Eq. (|4.6| ), taken e.g. at 
the right electrode, yields for the current of the bound state 

eA 2 



houn d {E) = e(|6| 2 - |/| 2 ) 



IE 



-D sin d>. 



(4.9) 



A useful formula for the current of the single bound state, which allows direct evaluation of 
the current from the bound state spectrum, is given by equation: 



14 



1(E) = 2e 



dE(4>) 



(4.10) 



where E((p) is the bound state energy band, Eq. ( |3.14| ). This formula is derived in Appendix 
B. Taking into account Eqs. 



the total current has the form [10-121531: 



eAD sin < 



: tanh • 



AJl — D sin 0/2 



(4.11) 



2y/l - L> sin 2 0/2 2T 
Thus, the dc Josephson current in tunnel junctions is carried only by the bound states, 
which is similar to the situation found in the other kinds of short weak links |p0ip4H26|j28|] . It 



follows from Eqs. ( [4.6|) , (|4.9| ) that the nonvanishing total current results from the imbalance 
of the bound state currents due to a difference in the equilibrium population numbers. 
Creation of a nonequilibrium population gives rise to a possibility to control the Josephson 
transport pTlMUl. 



V. INELASTIC SCATTERING 

Now we proceed to a discussion of inelastic scattering in voltage biased junctions. Ac- 
cording to our assumption $ = 0, explained in Sec. II, the applied voltage drop V is confined 
to the constriction; we also neglect, in order not to complicate the problem, a small time 
dependent voltage induced across the junction by the ac Josephson current (self-coupling 
effect ||54|| ). This implies the following dependence on time of the phase difference: 

(p = (po + 2eVt. (5.1) 

The appearance of factors with periodic time dependence in the boundary condition, Eqs. 
( p.8| ), (|2.9|) , gives rise to a more complex structure of the scattering states than in Eq. fl3.5|) . 
In order to satisfy the boundary condition, the outgoing part of the scattering states in 
Eq.(|3.5|) is to be constructed from the eigenstates of equation Eq. ( |2T7| ) with different 
energies E n = E — neV shifted with respect to the energy E of the incoming wave with an 
integer — oo < n < oo (sideband structure) 
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fe ) (0)= ( \u- E e.- iE * + V [ ] »* t -« (5.2a) 

(3 (o) =(!:)"^ E,+ ?(;),^ iE "' <5 ' 2b » 

For brevity we use the notation u n = ue„- While the incoming state is itinerant, the outgoing 
states can be either itinerant [Eq. ( |3.1| ) if \E n \ > A] or bound [Eq. ( |3.16| ) if \E n \ < A]. It is 



convenient to combine both the equations for the functions u n in a single analytical form: 

(5.3) 



/ „± 7n /2 \ 



\J2 cosh T n 



e 



A [ i(r nV /A 2 - El \E n \ < A. 

To find the scattering amplitudes in Eq. ( |5.2| ) we consider the boundary condition Eq. fl2.8| ) . 
It is important to mention, that this boundary condition was derived neglecting the energy 
dispersion of the normal electron scattering amplitudes d and r, which means that now this 
assumption should be valid for the entire interval of relevant energies E n . Let us first discuss 
j — 1 (hole-like quasiparticle coming from the left): 

1\ (a\ fc\ d /l + a z \/f 

0/ V 6 /, „ \f 2 V 1-aJ \c 

f\ 

I «n+l- ( 5 - 5 ) 

/ l,n+l 

It is convenient to separate the equations for normal and Andreev scattering amplitudes 



/ 1,71 


\/ / l,n 


d 1 


<\-o z 


+ 2( 






V l + a z 



l,n-l 



in Eq. (|5.5|) again using a procedure similar to Eq. (|3.8|) . The equation for the normal 
scattering amplitudes then becomes 

rc 1>n + (d/2)(V- n+1 f 1>n+1 + V+ n _J 1>n ^) = 5 n0 (5.6) 

rfl,n + (d/2){V+ n+l C hn+1 + Kn-lCl,n-l) = 0, 

16 



where the coefficients 



have the explicit form 



e-^m)/2 / cosh r n 
^ = ~ sinh 7n Vc^hTV (5 - 8a) 



m 



Km = cr n^m^-r \/ ' ( 5 - 8b ) 

sinn 7„ V cos h 1 m 

The equation for the Andreev scattering amplitudes reads 

a 1>n = (d/2)(C/-„ +1 / lin+1 + t/^.Ji.n-i) (5.9) 

where the coefficients are defined as 

= '"f^* 1 ^ , (5.io) 

(«+*, (7 2 u+) 

and have the explicit forms 



+ e (7n-7m)/2 / cos h Y n 

nm sinh 7 n V cosh T m ' ^ ^ 



^nm = - ff nffm r~r W 7"=-- (5.11b) 

sinh 7 n V cos h 1 m 



As can be seen from Eqs. (|5.6| ),( ^9|) , the inelastic scattering possesses a specific asymmetry: 
the forward scattered waves have odd side band indices and backward scattered waves have 
even side band indices, as illustrated in Fig. 4. Correspondingly, bound states with odd or 
even side band indices are induced either in the right or in the left electrode. We note that 
the scattering to any side band consists of both normal and Andreev components. 

It is instructive to compare the superconducting scattering diagram in Fig. 4 with 
the scattering diagram of normal junctions. In the normal limit A = 0, all the Andreev 
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amplitudes in Eq. Q5.9Q vanish [U^ = in Eq. (|5.11|) 1 and equations Eq. ( |5.6|) split due to 



V + = in Eq. (|5.8| ), which yields f n = c n _i = for all n ^ 1. Thus, the side band diagram 
in Fig. 4 reduces to the elementary fragment shown in Fig. 5a. This fragment corresponds 
to the scattering of a true hole, meaning a particle with spectrum Eh = —(p 2 /2m — jj) 
according to the BdG equation Eqs. (^.2|),(PT5|). In the ground state, T = 0, these holes fill 
all positive energy states E > while the negative energy states are empty. For the electrons, 
the corresponding diagram is sketched in Fig. 5b. In this diagram the chemical potentials 
in both electrodes are equal, while the energies of the incident and transmitted states are 
shifted by eV. This difference from the conventional diagram of normal electron tunneling 
in Fig. 5c (where the chemical potentials in the electrodes are shifted relative to each other, 
while the scattering is elastic) appears after separating out the superconducting phase in 
Eq. (p.4j) ; the conventional picture with shifted chemical potential can be restored by means 
of the gauge transformation of the normal electron wave function if) — > exp(— ieVt)if>. 

When superconductivity is switched on, A ^ 0, the incoming quasiparticle consists of 
both electron and hole components, and therefore the scattering diagram is a combination 
of the diagrams in Figs. 5a,b. One has also to take into account electron-hole conversions, 
which lead to the appearance of both electron and hole components in the upper as well 
as lower transmitted state. Continuation of this process creates the whole superconducting 
scattering diagram in Fig. 4. 

From a mathematical point of view, Eqs. (|5.6|) ,( [5l9D for the scattering amplitudes are 
second order difference equations which can not be solved exactly, except in some particular 
cases, e.g. a fully transparent constriction (r = 0) where equation Eq. (|5.6|) reduces to a 



binary relation p4| . In the general case, it is only possible to find asymptotical solutions 
using some small parameter. In the present case of a tunnel junction, there is a natural small 
parameter - the transparency of the tunnel barrier: Del. However, a straightforward 
perturbation expansion with respect to this parameter gives rise to divergences, similar to 
difficulties of the multiparticle tunneling theory (MPT) ]4^ , |57|J5^| ] . In order to formulate an 
improved perturbation procedure, it is convenient to rewrite Eqs. (|5.6p , (|5.9D in terms of the 



parameter A = D/4R, the true small parameter of the theory, as will be seen later. To this 
end we introduce new scattering amplitudes 

X k X k d* 

Cl,±2fc = — C±2k, /l,±(2Jfc+l) = ~2^~/±(2fc+l)) a l,±2k = X d± 2 k, (5-12) 

X k d 

6i,±(3fc+i) = ^r^±(2fc+i), 

which satisfy the equations: 

c « + XV~ n+1 f n+1 + V+ n _J n -i = (5.13a) 

a n = AC/^ +1 / n+ i + C^„_i/n-i (5.13b) 

6 n = AL r r j t n+1 c n+ i + L r „ n „ 1 c n _i, 

for n > 0. For n < one has to make the change V^ m —* V^ n _ m U^ m — > U^ n _ m , (n, m > 0) 
in the above equations. The equation for n = then reads 

c + \(V 01 f 1 + V + _ 1 f- 1 ) = l (5.13c) 

Let us now turn to the second scattering case in Eq. (|5.2|), j = 2 (hole-like quasiparticle 
incoming from the right). According to the symmetry relations of Eq. ( |3.6p the scattering 
amplitudes j = 2 differ from the scattering amplitudes j = 1 by cf) — > —<j), which in the 
present time-dependent case means that n ± 1 — > n =F 1- Taking into account this symmetry 
and also the property of scattering amplitudes in Eq. Q3.61 ), we introduce new scattering 
amplitudes 

X k \ k d* — 

fi,±ik = —c±2k, C2,±(2jfc+i) = /±(2fe+i)> ^2,±2fc = A a±2fc, (5-14) 

X k d- 

02,±(2Jfc+l) — ^70±(2fc+l), 

which satisfy the following equations (for n > 0): 
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C n + AKl+l/n+1 + L-l/n-1 = (5.15a) 

a n = XU^ n+1 f n+1 + U nn _ 1 f n _ 1 (5.15b) 

b n = \U nn+l c n+ i + ?7 nn _ 1 c n _i, 

Co + A(K,1/ 1 + V -i/-i) = 1. (5.15c) 

Equations (|5.15 ) differ from equations Eq. (|5.13| ) by 

V ± V T , U* -> f/ T . (5.16) 

In the case of electron-like quasiparticles incoming from the left, j = 3, the symmetry of 
Eq. (|3.7| ) involves transformation 7 — > —7, which means transformation of the coefficients 
— >■ —(Jn a -mY^ ■, Z7 ± — > —a n a m U T in Eqs. ( |5.6|) , ( |5.9| ). This allows us to relate the 
scattering amplitudes of this case to the solutions of Eqs. ( |5.15| ): 
X k X k d — 

«3,±2fc = — &±2kC±2k, &3,±(2fc+l) = — ^^" Cr ±(2fc+1)/±(2A:+1)5 c 3,±2fc = A a ±2 k^±2k, (5-17) 

/s,±(2fc+l) = ^-^^±(2k+l)b±(2k+l)- 

In a similar way the scattering amplitudes of electron-like quasiparticles incoming from the 
right, j = 4, are related to the solutions of Eqs. ( |5.13| ). 

,(4) A^ (4) p f U) _ \k (r -iq\ 

b ±2k — — a ±2kC±2k, a ±(2k+l) — ~-^ a ±(2k+l)J±(2k+l), J±2k ~ A &±2ka±2k, (0.18J 



C ±(2fc+1) ^-r^±(2fc+l)fr±(2fc+l). 

According to a symmetry of the coefficients in Eqs. ( 5.13|) , ( |5.15| ). 



X k d* 



V± m (-E) = V±:_ m (E), U± m (-E) = U±* n _ m (E), (5.19) 

all scattering amplitudes with positive and negative incoming energies are related through 

a n (-E)=a*_ n (E), (5.20) 
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and similarly for the other amplitudes. 

Let us now present a formal solution of Eq. ( |5.13| ) for n > on the following form |55 



2k+l 

fa+i = (-1)" II ^c , (5.21) 

1=0 



where the quantities Si are defined as 



<~>2fc — ~ 7 ) >->2fc+l — , [p. II) 

J2k-1 c 2k 



and satisfy the recurrence relations 



V + V~ 

c v 2k,2k-l Q v 2k+l,2k ( r qq\ 

1 + AV/ 2fc,2A:+1^2fc+l 1 + A * / 2fc+l,2A : +2^2fc+2 



The quantity Co in Eq. (|5.27Q is given by 



c °- i+Kv^+vi^y (5 - 24) 



It is convenient to express the functions S n in Eq. ( |5.23| ) through the relation S n 
V^n+\/Z n , were the denominators Z n (n ^ 0) satisfy the recurrence 



^n = l + A^±i, a±=4r— , (5.25) 
Z n+1 smh 7 n 



(± corresponds to even/odd n), and to define Z as the denominator of c , Eq. (|5.24j) : 

Z =l + + A^I. (5.26) 

Z\ Z_i 

Using the above notation, one is able to express the coefficients of the normal forward 
scattering, |/ n | 2 , on the form 

e 70 n i 

cosh 70 1 Z 1^ | Zi smh 7; | 2 

provided a~ 7^ 0. The equation for the coefficients of the normal backward scattering, \c n \ 2 , 
differs from Eq. ( |5.27| ) by e r " — > e~ Vn . The relation between the amplitudes of the Andreev 
and normal forward scattering in Eq. ( |5.13| ), taking into account Eqs.( |5.22"| ), ( |5.23| ), Q5.25| ), 



has the form 
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7n fl - A- 



2e -7n + l 



fn- 



sinh7 n+ iZ„ + i 

In a similar way, one can express the solution of Eq. (|5.15 ) for n > on the form 



(5.28) 



l/n| 



-70 



cosh7o|Z | 



re r "coshr n ]J 



^ |Z;sinh7i 



b n = -e 7n (l + A— 

\ SIT 



|2 ' 



sinh7 n+ iZ n+ i 



fn: 

(5.29) 



where 



1 + A- 



aja 



n u n+l 



1 + A- 



a a x 



+ A- 



"-1 



(5.30) 



We note that Eqs. ( |5.29| ), ( |5.30| ) differ from Eqs. ( |5.27| ), ( |5.28| ) by 7 n — >■ — 7 n everywhere. 

Equations for the scattering amplitudes with negative side band indices, n < 0, can be 
derived in a similar way, and the result differs from the above equations for positive side 
band indices Eqs. ( |5.25| )-( |5T3"0| ) by the substitution 



In -> — T-Hj n ^ 0, 



(5.31) 



introduced everywhere except in Z Q and Z . 



VI. QUASIPARTICLE CURRENT 

In the nonstationary problem under consideration, the density matrix determining the 
current Eq. (|4.1|) is time dependent, and its dynamic evolution can be described by an 
equation similar to Eq. (|4.3|), 

(*(f,t) ¥^,t))=J2^x(r,t)f x ¥ x ^,t), (6.1) 

A 

ty\ are now solutions of the time dependent problem, Eq. ( |2.5| ), with initial conditions 
corresponding to the eigenstates of the initial Hamiltonian with the eigenvalues A, and 
occupation numbers fx of these initial states. We consider inelastic scattering states, Eqs. 
([2.6fl, (|5.2|) , as the propagators ^\(t) in Eq. ( |6TT| ) with A corresponding to the complete set of 
the incoming states A = (E,j); according to the assumption about local equilibrium within 
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the electrodes, the incoming states possess the Fermi distribution of occupation numbers, 
/e? = n F (—E). Thus the current Eq. ( |4.1| ) takes the form: 



r rl h' I h'\ 00 00 

m = -ej J^rM-E) E j NeVt E Y,P(^(E,n),^(E,N + n)). (6.2) 



dE\E\ 

!\E\>A 2tt£ _ v __ x „__ x ; , 



The current in Eq. ( |6.2j ) consists of a time independent part, N = 0, which is formed 
by incoherent contributions of all the side bands (the quasiparticle current), and of a time 
dependent part, N ^ which results from interference among the different side bands (ac 
Josephson current). The difference between the side band indices N is an even number 
since the side band index is either even or odd depending on the electrode; therefore the 
time-dependent current oscillates with the Josephson frequency uj = 2eV. 

In this paper we will concentrate on an analysis of the time-independent quasiparticle 
current. Similarly to Eq. (|4.7f ) we will calculate the current using the transmitted states, 

/ = ~ / , A d -^n F {-E) £ [ E (l/inl 2 - IM 2 ) + E (kni 2 - M 2 )} . (6.3) 

Z7YJ\E\>A n=-oo L i=l,3 j=2,4 J 

Using the scattering amplitudes introduced in the previous section through Eqs. ( |5.12|) , 
( |5.14j ), ( ^.17|) and ( |5.18|) we express the current in Eq. (|6.3|) in the following form: 



/ = - / , A ~~}~~ n F{—E) E(^n - K n ), (6.4) 

where 

K n = \\ n \R- l \f n \ 2 - \b n \ 2 ), K n = X^iR'^JJ 2 - \b n \ 2 ) = K n (-j). (6.5) 

The factor of two appears in Eq. ( p.4| ) due to equality of currents I\ and J4 and of currents 
I2 and ^3 in Eq. (|4.8|) , equalities which hold also in the nonstationary case. However, there 
is no balance between currents of these two pairs any more. The symmetry of Eq. (5.20) 
allows us to reduce the interval of integration in Eq. ( |6.4j ) to the semiaxis E > 0, 

I=~ r^tanh^E(^-^n). (6.6) 

n JA $ 2T odd 

The side band currents K n in Eq. ( Jt>.5|j are proportional to the powers of the small parameter 
A, K n ~ A™. Therefore Eqs. ( |6.6|) , (|6.5|) present a perturbative expansion of the current, 
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convenient for analysis in the limit of low barrier transparency. The following sections are 
devoted to such an analysis of the structure of the current in Eq. ( |6.6| ) . 



VII. EXCESS CURRENT AT LARGE BIAS 

To make some observations useful for analysis of the subgap current, it is instructive 
first to discuss the simpler case of large bias eV ^> A, which is well studied in literature 
17| , pT|p3| , p!2[] . Simultaneously we will derive the explicit analytical expression for the current 



in this limit valid in a whole range of the junction transparency < D < 1. The asymptotic 



expansion of the current with respect to the small parameter A/eV has the form [|L7] : 



e 2 DV ( A 

1 = + I exc (D) + (7.1) 



7r ' \eV 

where the first term is the tunnel current of the normal junction and the second term is a 
voltage-independent excess current which represents the leading superconducting correction. 

A main simplification in this case is that the side band currents K n and K n , \n\ > 1 
diminish when the bias voltage increases. This follows from an estimate of the transmission 
amplitudes in Eqs. (|5.27|) - (|5.29|) , which contain products of factors | sinh ^ | — 2 which are 



small at large voltages, |sinli7 fc |~ 2 ~ (A/eV) 2 , because of the large interval of involved 
energies, E ~ eV. Furthermore, inspection of the amplitudes /_i and fi shows that they 
are also small due to the factors e~ 7o ~ 71 ; therefore the non-vanishing part of the current Eq. 
( |6.6|) in the limit eV ^> A becomes 

e r°° dE E E — . . 

/ = -/ tanh— (7.2) 

The essential fragments of the scattering diagram in the large bias limit are shown in Fig. 
6. 

The structure of the current in Eq. ( |7.2| ) is essentially determined by the presence of a 
gap in the spectrum of the side band n = 1; this causes different analytical forms of the 
current Ki in the regions \E\ < A and \E\ > A. We note that the spectrum of the side band 
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n = — 1 possesses no gap: E-% > A for E > A. According to this we divide the integral in 
Eq. ( |7.2| ) into three parts: 



/< + I A + I 



>■ 



The first part corresponds to the current of the states in the side band n = 1, which lie 
below the gap, E\ < —A. The second part corresponds to the current of the states of the 
same side band lying within the gap, —A < E\ < A. The third part combines contributions 
from the remaining states of the side band n = 1, A < E±, and from the all states of the 
side band n — — 1. Making use of the approximations 

\Z \ 2 = \Z \ 2 « |l + 2A(£ + 0/£| 2 , \Z-i\ 2 = \Z X \ 2 = \Z 2 \ 2 = |Z_ 2 | 2 « 1, (7.3) 

\Z X \ 2 = \Z^\ 2 = \Z_ 2 \ 2 = \Z 2 \ 2 « l/R 2 , 
it is possible to express the integral 7< on the form (we restrict ourselves to the limit T = 0): 



n Ja £ nR Ja iZl ttR J a C^o \ £ / 

where the limit of integration in the last term is extended to infinity since the main contri- 
bution to this integral comes from the energies E ~ A <C eV. Separating out the normal 
junction current we may express Eq. (|7.4|) in the following form 



e 2 DV 8e\ 2 r°° dE 

I< - 



EA 2 , , E(E — £) 1 



(7.5) 



7T 7T 7a Z q 2 £ 

We note that this current is always smaller than the normal current. It is convenient to 
express the integral I a as 

e r ev+A 7? 16eA 2 r A A 2 
Ia = ~ dE^K 1 = ^ dEj^-, (7.6) 

71 JeV-A £ 71 JO \<iZo\ 

obtained by using the relations between the functions Z n [which result from their definition 
in Eq. ( g^g) ] 

Z n (E + eV) = Z n -!(E), Z (E + eV)Z 1 (E + eV) = Z (E)Z^(E). (7.7) 
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Inspection of the equation for / 



p poo pi e poo pi 

/> = - / dE-K x - - / dE-K- U (7.8) 

71 JeV+A 4 71 J A £ 

shows that both integrals diverge at the upper limit E = oo, which means that the states 
lying far from the Fermi level formally contribute to the current, while the quasiclassical 
approximation of Eq. ( |2.6|) assumes, that all relevant states lie close to the Fermi level. To 
get rid of this formal divergence one commonly shifts by eV the variable in the first integral 
in Eq. ( |7.8|) . Using again the relations (|7.7|) we may express this integral on the form 



where the first term has the same analytical form but the opposite sign compared to the 
divergent term in the second integral in Eq. (|7.8|) , 

« r dE ^,^r dE (E-o 



7i R J a tZoZ-t nR J a iZl 
After elimination of the divergent terms, the integral in Eq. ([7.8Q takes the form 



The currents in Eq. ( |7.9|) and Eq. ( |7.6|) are positive and overcompensate the missing part 
of the current in Eq. ( |7.5|) . Collecting Eqs. (|7.5| ), (|7.6| ) and (|7.9| ) we find after some algebra 
the following explicit equation for the excess current in Eq. (|7.1|), 



\<oek\ 2 R 



7Y 



D 2 1 + VR 

1 = In 



(7.10) 



2(1 + r)Vr l-VR}' 
which is valid in the whole interval of junction transparency < D < 1. Asymptotics of 



this expression coincide with the results presented in literature [|TT| , p3f both in the limit of 
fully transparent (D = 1) constrictions, I exc = 8eA/37r, and in the limit of low-transparency 
(D <^ 1) tunnel junctions, I exc = eAD 2 /7i. 

The above calculation reveals an important difference between the structure of the cur- 
rent in normal and superconducting junctions. In normal junctions, the current, e.g. in the 
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right electrode, see Fig. 5c, results from scattering states lying above the local chemical 
potential, E > fi — eV, while contribution from the energy interval E < fi — eV is equal to 
zero due to mutual cancellation of currents of the scattering states incident from the left and 
from the right (in Fig. 5a the current carrying energy region corresponds to negative energies 
Eh < 0). Thus the total current coincides with the current of real excitations emitted from 
the contact, which is consistent with the nonequilibrium origin of the current in the voltage 
biased junctions. In superconducting junctions, only " across-the-gap" current J< is clearly 
related to the real excitations emitted at the right side of the junction where the current 
is calculated (Fig. 6a) - the dissipative character of the currents /a and J> is not obvious. 
However, one should take into account the creation of real excitations at the left side of the 
junction via back scattering into the side band n = 2 (Fig. 6b, c). Although the current of 
this side band exists only at the left side of the scattering diagram it should have an effect 



at the right side due to continuity of the current at the interface, Eq. ( f4.5|) , and therefore 
it should be distributed among the states of the side band n — 1. As our calculations show, 
this "kick" current partially flows through the Andreev bound states, involving the current 
Ja, Fig. 6b, which convert this current into a supercurrent outside the junction. It is also 
partially distributed among the scattering states with positive energies, current J>, Fig. 6c, 
in the form of imbalanced ground state currents. 

VIII. SUBGAP CURRENT 

In this section we turn to a discussion of the tunnel current in the subgap region eV < 2 A. 
A basic property of the subgap current is the presence of temperature independent structures 
on I-V characteristics - the subharmonic gap structure (SGS). SGS in tunnel junctions was 
discovered in experiments by Taylor and Burstein [pUfl , and the first theoretical explanation 
was given by Schrieffer and Wilkins in terms of multiparticle tunneling (MPT) . Recently 



SGS has been observed in a number of experiments on transmissive tunnel junctions 
Although there is a possibility to attribute SGS in planar junctions to normal shorts, the 
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observation of SGS in superconducting controllable break junctions || provided convincing 
confirmation of the existence of SGS in the true tunnel regime. 

The existence of SGS in tunnel current can be established within the MPT theory 
by means of rather simple perturbative arguments P2| , |57| -p9|1 . On the basis of the tun- 
nel Hamiltonian model, assuming a small perturbative coupling between electrodes, one 
can calculate the probability of tunneling in n-th order of perturbation theory. Such a 
probability is proportional to a product of filling factors of the initial and the final states: 
Uf(E)[1 — np{E — neV)}. At zero temperature this factor is equal to zero outside the inter- 
val A < E < neV — A, which selects the quasiparticle transitions across the gap, i.e. the 
processes of creation of real excitations relevant for the tunnel current. Such a restriction 
places the threshold of the n-th order current at eV = 2A/n, and a sequence of current 
onsets of magnitudes ~ D n at the voltages eV = 2A/n forms the SGS of the tunnel current 

IIS- 

In our approach, the filling factors of final states do not enter the equation for the current 
Eq. fl6.6|), and the existence of SGS is therefore not obvious, although the side band currents 
Eq. ( |6.5D successively diminish with increasing side band index. However, attribution of 
the nonequilibrium tunnel current in biased junctions to the current of real excitations is a 
general physical argument which should be automatically met in any correct theory. In fact, 
and this also follows from the discussion of the previous section, the true tunnel current is 
hidden in Eq. ( |6.6|) : it results from partial cancellation of large contributions of different 
scattering modes. The cancellation is nontrivial because of mixture of currents of different 
side bands, the odd side bands containing information about the currents of the even side 
bands and vice versa. This means that a finite perturbation expansion of Eq. ( |6.6| ) is not 
satisfactory and will not adequately correspond to the perturbative structure of the true 
tunnel current. To reveal such a structure one must rearrange the series in Eq. ( |6.6| ). 

To this end we consider a general term K n , n > in Eq. ( |6.6| ). It follows immediately 
from the explicit form of the normal and Andreev transmission coefficients, Eqs. ( |5.27| ) and 
( |5.28| ), that the leading term with respect to A in K n is proportional to a factor (1 — e 2Tn ) 
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which is equal to zero if \E n \ < A. Having made this observation we express the quantity 
K n in the following form: 



n+l -2IY, 



R 



\ n 6(E 2 n - A 2 )e^ sinh 7n |/ n |^ + 4A n+i e 



fn 



J n+l 



n+l, 



n+l 



'n+l | 



+ R8 (f__£s±l)_A 



"7n+l 



sinh 7 



n+l 



(8.1a) 



(8.1b) 



sinh 7„+i 

In Eq. ( |S-laj ) the first term represents the main contribution of the n-th side band to the 
current: it is proportional to the probability of normal scattering to the n-th side band and 
it does not contain the contribution of the side band states lying inside the gap \E n \ < A. 
Using the recurrence relation (|5.25 ) and remembering that A = D/AR, after some algebra 
the function F n in Eq. ( |3.1b| ) becomes 



-6(E 2 n -A 2 )—^— 
R K n ; tanh 7r , 



A 



sinli7 n Z n+ i 



Gn 



+1, 



G 



n+l 



J n+l| 



M ■ , + ) - A 



oln+l 



sinh j. 



n+l 



(8.2a) 



(8.2b) 



sinh 7 n+ i 

Substituting Eq. ( p.2|) into Eq. (|8.1a| ), we find that the second term in the equation for K n , 
proportional to A n+1 , has analytical structure similar to the first term in the same equation, 
proportional to A n , namely, it consists of the probability of normal scattering to the (n+ 1)- 
th side band [cf. Eq. ( |5.22|) 1 and it does not include the contribution of the side band states 
lying inside the gap |J5 n +i| < This allows us to associate this term with the effective 
contribution of the nearest even side band. 

A similar transformation of the function G n+ \ in Eq. (|8.2|) yields the recurrence: 



G 



n+l 



^ ( ^- A2) t^7^ 



A 



n+2 



(8.3) 



sinh 7^+1^+2 

Combination of Eqs. (|8.1a|) - (|8.3|) shows that the next term of the current K n , proportional 
to A n+2 , has the same analytical structure as the leading term in the current K n+2 of the 
next odd side band, and therefore it can be regarded as a renormalization of that current. 
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Continuing this procedure by systematic use of the recurrence Eqs. 
the following expansion for the current K n in Eq. ( |6.5| ): 



2A n 



Q{El - A 2 )Q n 



4A 



n+i 



6{E 2 +l - A 2 )e~ r " coshr n g n+1 



R 



~M) we obtain 



(8.4) 



4A 



n+2 



R 



-6(E 2 n+2 - A 2 )e 



^ r ri -i-2r\ 1 . 



- 1 cosh T n Q n+2 



4A 



n+3 



R 



-9(E 2 +3 - A 2 )e 



2\ r„+2r n _|_i— 2r„^ 



cosh V n Q n+3 + ... 



where we have introduced the quantity Q n defined for all n as 

1 



,70 



Q, 



■ sinh 7„ cosh T n ]J[ 



■5) 



cosh7o|Z | 2 " ^|Z z sinh7;| 2 

One can derive similar expansions for the currents K n as well as for the currents of the 
side bands with negative n < 0. Expanding each term of the series in Eq. (|6.6|) using Eq. 
( p.4[) and collecting terms with the same factor A™, we can finally express the series in the 
form 



odd nj^O 



(8.6) 



The last summation is done over all odd and even integer n, and the renormalized coefficients 
have the form 



K n = \ n 6(E 2 - A 2 )(AQ n /R) (1/2) + coshl^e" 1 " 



-r„,_o+2r 



2-I-Z1 n _i 



(8.7a) 



+ coshr„_ 4 e- r "- 4+2r "- 3 - 2r "- 2+2r "- 1 + ... + coshl^e"^^-^^-^' 1 - 1 



for odd n > 0, and 

K n = X n 0{E 2 n - A 2 ){AQ n /R) [coshlVie-^- 1 + coshr^-^- 3 ^"- 2 " 2 ^- 1 + ... (8.7b) 

+ coshr 1 e- ri+2r2 - 2r3+ -- 2r "- 1 " 

for even n > 0. 

The representation of Eqs. ( |8.6| ),( |8T7| ) is exact. One can regard a general term of the 



series as an effective renormalized current of the n-th side band. In fact, this effective current 
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consists of the contributions of all side bands with odd indices smaller then n. An important 
feature of this representation is the presence of the ^-function in the general term, which 
allows us to separate out in Eq. (|6.6|) the part of the current which is obviously responsible 
for the SGS, 

00 g pneV-A J? E ~ 

Isgs = J2- dE - tanh— (K n -K n ) (8.8) 

n=l n Ja S 2T 

One might expect (cf. Ref. |55| ) that Eq. (|8.8| ) represents the subgap tunnel current at zero 
temperature and that the remaining part of the current of Eq. ( |6.6| ), 

00 e 



h = I ~ I SGS = E - 



POO f-f] j-jj POO j^J f-jj ~ 

/ dE - tanh —(K n -K n )+ dE - tanh — {K_ n - K_ n ) 

(8.9) 



corresponds to the current of thermal excitations. However, this separation is not exact: 
an analysis shows that the current in Eq. (|8.9|) does not vanish completely at T = 0, but 
contributes a small residual part. An important property of this residual current is that 
it does not contain any structureless component but demonstrates behaviour similar to the 
current Isgs in Eq- (|8.8|), thus resulting in a small correction to Eq. 



IX. SUBHARMONIC GAP STRUCTURE. 

The explicit analytical expressions ( |8.8| ) and ( |8.9| ) provide a basis for numerical calcu- 
lation of the subgap current for small A (low transparency) with any desirable accuracy. 
However, they are also convenient for qualitative discussion of the SGS. In this section we 
will analyze the SGS at zero temperature on the basis of Eq. ( jSJgp . 

The current-voltage characteristic Isgs{V) in E<1- ( |§.8|) has a complex form consisting 
of a sum of renormalized side band currents I n (V, A): 

e rneV-A E ~ — 

Ml - A). Ml .A) '■• • 

n=l 

The partial current- voltage characteristics I n (V,X) are similar to each other, and it is con- 
venient to analyze them independently. 
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00 g rneV-A E ~ — 

Isgs(V, A) = E W A )> W A) = - / dE -(K n - K n ). (9.1) 

ZZ^ TV J A t 



According to Eq. ( |9.1| ) the partial current I n starts with an onset at the threshold voltage 
V n = 2 A/en. In the limit A — > the onset is infinitely sharp and the magnitude of the onset 
is 

2ti 71 2n 

/ n (K,A-0) = eA^— — . (9.2) 

The jumps of the current at the thresholds result from the singular denominators in Eqs. 
Q5.271 ), (15.29| ), related to the singular density of states at the side band energy gap edges, 
sinli7fc = 0. Accumulation of these singularities in the high order scattering amplitudes 
yields a tremendous increase of the partial currents well above the corresponding thresholds 
- this is what causes the failure of multiparticle tunneling theory f57H59"l . In our theory, the 
singularities are regularized by the factors 

n 

Pn=U \Z k \ 2 (9.3) 

fc=0 



in the denominators of the scattering amplitudes Eq. (|5.27 ). These factors are expressed 
through the continued fractions Z n , Eq. (|5.25| ), which therefore should be calculated with 
sufficient accuracy to preserve the singular parts of Z n which provide regularization of the 
integrals. 

The first order current I\ in Eq. ( |9.1| ) corresponds to direct one particle scattering to 
the side band n = 1, Fig 7a. The explicit form of the current I\ is 

Il = *4r v - A dE m(E±i + E^), (9 . 4) 

irflii «, \ -Pi Pi J ^ ' 

In the limit A — > this current coincides with the quasiparticle current of the tunnel Hamil- 



tonian model pll , |54]| . At finite A the threshold onset of the current at V = V\ is washed 
out. To evaluate the width of the onset we truncate the continued fraction in Pi assuming 
Z-i = Z2 = 1, obtaining 

Pi w |(1 + Ado cQXl + Xa^a^) + Xa^a^\ 2 . (9.5) 

The function Pi has a similar form. The regularization effect of the threshold singularity is 
provided by the most singular term Xa^af in Eq. ( |9.5p . Keeping this term we obtain in the 
vicinity of the threshold, e(V — V\) A, the following result 
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t n/\ 2eAA / eV-eVA ^ Aa sin 2 9 

h{v) = ~^R f l^^J ' f{z) = Jo de PTW (9 - 6) 

According to this formula the onset width is e(V — V\) ~ AA. 

The second order current I 2 corresponds to the creation of a real excitation during 
quasiparticle backscattering into the side band n = 2, Fig. 7b, and appears as the current of 
transmitted states of the side band n = 1 (cf. the excess current in Sec. VII). In the vicinity 
of the threshold, V 2 < V < Vi, this current exists only in the form of currents through 
the bound states and therefore it is completely converted into supercurrent far away from 
the junction. At larger voltages, V > V\, the side band n — 1 extends outside the energy 
gap (see Fig. 8a), which also makes the current I 2 partially consist of contributions from 
extended states. The explicit expression for the second order current is 

4eA 3 A 2 r 2eV ~ A Irn \E 2 \ , „ /e" 70+ri e^- ri \ , wN 

dE ' coshrJ - +^=— . (9.7) 



tiR Ja e^Kil 2 \ P2 Pi 

Omitting the A-dependence of P 2 in Eq. ( |9.7| ), one gets the two-particle tunnel current of 
Schrieffer and Wilkins p2]j57[| . To keep the singular terms in P 2 one has to truncate the 
continued fractions in Eq. ( |5.25| ) assuming Z_i = Z% = 1, which yields 



P 2 w |(1 + AcQao )(1 + Xa^a 2 ) + Aa^a^(l + Aa^a^")| 2 . (9.8) 



The threshold singularity results from the small product ££2 hi denominator of Eq. 
However, there are no singular terms in Eq. (|9.8j ) proportional to aoa 2 among the terms 
linear in A. Such terms are quadratic in A and they provide, along with the terms Aao and 
Xa 2 , the width of the onset: e(V — V 2 ) ~ A 2 A. This onset is sharper than the onset of the 
current I\. 

The threshold singularity in the current I 2 is typical for all higher order currents n > 1. 

The appearance of the first side band outside the energy gap at V = V\ is manifested 
through a spike in the current I 2 . Indeed, if V « V\, the nodes of £1 overlap the nodes of £ 
and £2 at the lower {E = A) and at the upper (E = 3A) limits of integration in Eq. ( |9.7|) 
respectively (see Fig. 8a). This singularity yields an increase of the current I 2 when the 
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voltage approaches Vi, 



I 2 ~ eAA 2 



A 



eW-vy 



Regularization of the integral, which is provided by the singular terms ActiOo and \a1a2 in 
Eq. (|9.8|) at the lower the and upper integration limits respectively, yields 



Further analysis shows that maximal magnitude of the current is achieved slightly above 
V = Vi, after which the current rapidly decreases (see Fig. 9). At voltages V > V\ the 
singular point £1 = remains within the integration region, which gives rise to enhancement 
of the magnitude of the current by a logarithmic factor in comparison with the current 
magnitude near the threshold V2, 



At large voltage V ^> V\ the current I2 forms the excess current Eq. (|7.10|) . It is interesting 
to note, that in this limit the logarithmic factor is compensated for by the current I r , Eq. 
QB.9|) which yields the A 2 -dependence of the excess current. 

The third order current -Z3 at voltages close to the threshold V3 results from the combi- 
nation of one-particle tunneling into the side band n = 3 and excitation of the transmitted 
Andreev bound states of the side band n = 1, Fig. 7c. The probabilities of these two pro- 
cesses are related as 1:2 at threshold, Eq. ( |8.7| ). Successive emergence of the bound states 
of the side bands n = 1 and n = 2 outside the energy gap at V = V2 and V = Vi, Fig. 8b, 
gives rise to the current peaks. The current ^3 has the explicit form 



h(v 2 ) ~ vT 



h(V > Vi) 



eAA 2 In A 
R 



(9.9) 




-ri+2r 2 



) e- 70 (l + 2coshrie ri - 2r2 ) 



^3 J 
(9.10) 



with the regularization factor 



P 3 « |(1 + Aa_ 1 a )(1 + \a 1 a 2 )(1 + Aa 3 a 4 ) + AaJ" af(l + Xa^a^)] 2 . 



(9.11) 
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The current peak at V = V 2 results from the overlap of nodes of £ and £2 at E = A and nodes 
of £1 and £3 at E = 2A, similarly to the peak of the current I 2 - These singularities yield 
again an increase of the current inversely proportional to the square root of the departure 



from the voltage V 2 : I3 ~ eA 3 ^ 2 X 3 / y e(V 2 — V). However, since the factor P 3 , Eq. ( p. lip , 
contains neither a term \a$a2 nor a term \a1a3, regularization of the singularity is provided, 
e.g. at E = A, by the terms Xa or Aa 2 , which gives rise to a more pronounced peak with 
magnitude 

h{V 2 ) 1 



/3(H) y (9 ' 12) 

We note, that the magnitude of this peak is comparable with the magnitude of onset of 
the current I 2 . The second peak at V — V\ results from the overlap of nodes of £2 and £3 
at E = 3 A; this gives rise to increase of the current J3 near the voltage V = V% inversely 
proportional to the first power of the distance to this voltage: J 3 ~ A 3 A 2 /(Vi - V). The 



divergence is regularized by the term Aa2CZ3 in Eq. ( |9.11| ) which results in a peak with 
magnitude 

h{Vi) 1 



(9.13) 



h(V 3 ) A' 

Thus the heights of both peaks of the current J3 are of the same order in A, although the 
peak at V « V\ is sharper. 

In a similar way, all of the high-order currents result, in the vicinity of their thresholds, 
either from Andreev bound state currents (even n) or from a combination of Andreev bound 
state currents and the current of a single real excitation (odd n). The number of excited 
Andreev states is correspondingly n/2 ot (n — 1)/2. Singularities similar to the singularity of 
the current J 3 at the voltage V « V% exist in all high order currents, where they cause even 
more pronounced current peaks because of the absence of terms Adfcdfc+i in the corresponding 
smearing functions P n . Due to this property, the magnitudes of such peaks exceed the 
threshold magnitude of the corresponding current by two orders of A: (I„)max ~ eA\ n ~ 2 /R. 

The above discussion reveals the current peaks to be essential features of the SGS of 
tunnel current in addition to the current onsets, Fig. 9 (these peaks are seen also in the 

35 



numerical results of Refs. |34|]3~9|1 ). It allows us to establish a general classification of singu- 
larities causing peaks in partial currents I n . They result from the overlap of singularities of 
the side band density of states. It is easy to see that the singularities of only two side bands 
can overlap, the condition of the overlap for m-th and /c-th side bands having the form 

E-keV = A, E - meV = -A. (9.14) 

This condition is met at voltages eV = 2A/(m — k) for all integer < k < m < n. 
The magnitude of the current peaks depends on whether the overlapping side bands are 
neighbours or not, and whether the side band index is inside or at the edge of the interval 
(0,n). 

I. m — k = 1, m = n or k = 0: edge-type singularity, neighbour side bands. This type 
of singularity forms the peak of the current I2 at the main threshold V\. The magnitude of 
the current peak is: (^Jmax ~ eAVX/R. 

II. m — k > 1, m = n or k = 0: edge-type singularity, non-neighbour side bands. This 
type of singularity forms the first peak of each current I n , n > 2 at voltage V n -\. The 
magnitude of the current peak is: (in) max ~ eAA n_1 /i?. 

III. m — k — 1, m < n, k > 0: internal singularity, neighbour side bands. This type of 
singularity forms the last peak of each current I n , n > 2 at voltage V\. The magnitude of 
the current peak is: (i n )max ~ eAX n ~ l /R. 

IV. m — k > 1, m < n, k > 0: internal singularity, non-neighbour side bands. This type 
of singularity forms all intermediate peaks of each current I n , n > 3. The magnitude of the 
current peaks are: (i n )max ~ eAX n ~ 2 /R. 

X. CONCLUSION 

In this paper we have considered superconductive tunneling as a scattering problem 
within the framework of Bogoliubov-de Gennes (BdG) quantum mechanics. An essential 
aspect of this problem is that the scatterer consists not only of the potential of the tunnel 
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barrier but also of the discontinuity of the phase of the order parameter. In equilibrium 
(zero bias, dc Josephson current) the scattering problem is elastic. The peculiar feature 
of the elastic scattering problem in short junctions, considered here, is that the balance of 
currents of the scattering modes is not violated: the supercurrent flows only through the 
superconducting bound states (for a more general discussion, see In the presence 



of voltage bias the scattering is inelastic due to the time dependence of the component 
of the scatterer related to the superconducting phase difference. Generally, the currents 
of all inelastic channels constitute together the components of the tunnel current flowing 
through the biased junction: the quasiparticle current corresponds to the incoherent part 
of the inelastic side band contributions, and the ac Josephson current corresponds to the 
interference of the side band contributions. 

There are three distinct components of the quasiparticle tunnel current at zero tempera- 
ture: (i) the current of quasiparticles excited above the ground state, (ii) the current trough 
Andreev bound states converted to supercurrent outside the junction, and (iii) the imbalance 
current of the ground state modes. At large bias voltage, eV ^> 2A, the first component 
corresponds to the single particle current of the normal junction, while the other compo- 
nents cause excess current. When voltage is decreased, redistribution of current among the 
components gives rise to subharmonic gap structure (SGS) in the form of current onsets 
and current peaks. Within the voltage intervals 2A/n < eV < 2A/(n — 1) with even n, the 
tunnel current entirely consists of currents through the Andreev bound states [component 
(ii); e.g. Fig. 7b], the states of all side bands with odd indices smaller than n contributing 
to the current. If n is odd, a real excitation current of the side band n [component (i); e.g. 
Fig. 7a,c] is also present in the tunnel current. Opening of new channels of tunneling of 
real excitations gives rise to current structures. Thus, SGS reveals the discrete nature of the 
side band spectrum. The structure is the more pronounced the smaller is the transparency 
of the junction. 

Since each Andreev state provides transfer of one Cooper pair through the junction for 
every incident quasiparticle, n particles will tunnel in the interval 2A/n < eV < 2A/(n — 1). 
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One may be surprised by the participation of a large number of bound Andreev states in 
current transport at low voltages: this appears to contradict the fact that subgap current 
diminishes at zero bias. After all, the probability of the scattering into side bands does not 
depend on the bias and is proportional to powers of D. The solution to this paradox results 
from increasing compensation of currents between the normal and Andreev channels in each 
side band with decreasing voltage, which gives the required voltage dependence of the total 
current. 
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APPENDIX A: BOUNDARY CONDITIONS 



The quasiclassical boundary condition in Eqs. (|2.8| ), ( |2.9| ) has been derived in Ref. |3T| 



using the method of Ref. ||21|| . Here we present simple arguments which lead to this boundary 
condition. We consider the more general case of asymmetric junction, using an asymmetric 
version of the Hamiltonian of Eq. (|2.2|) with the same restriction on the length of the 
non-superconducting region: L <C £ - We include a contact potential difference into the 
potential U(x) which implies that this potential may have non-vanishing asymptotic values 
at infinity: U(— oo) 7^ U(oo) 7^ 0. If the junction has more than one transverse transport 
mode we assume that these modes are not mixed. 

A one-dimensional quasiclassical wave function of a given transverse channel in the right 
electrode has the form [Eq. (|2.6|)1, 

ttflOM) = E -^e l Pf p * dx e^ 2 ^ R (x,t), (Al) 



p s/ v R 



38 



with similar expression for the left electrode. ip R are slowly varying two-component wave 
functions on the scale of 1/pr, where Pr(x) = y / 2m^(/i — Ur — E±r(x)). This equation is 
valid over the distance x ^> 1/pr from the junction, and in the spatial region 1/pr Ci <(o 
the functions ip^ are almost constant. 

From another point of view, at large distance from the junction, \x\ ^> 1/pr^l, the 
function \I/ can be expressed in the form of a linear combination of the scattering states at 
the Fermi level, 

= Cixi + C2X2 (A2) 



Xi 



{1/^/vRjde 



(A3a) 



x > 0, 



X2 



{l/y^)[e- ipLX + re ipR % x>0 



(A3b) 



'l/^JZ)de- ipLX , x<0. 
Comparing Eqs. (|A2|) , (|A3|) with Eq. (|A1|) in the region 1/ pp <C \x\ <C £0 we have 



d = e i<TzXL/2 ^, C 2 = e WzXfl/ V K) rCi + dC 2 = e ia * XL/2 ?p 



(A4) 



which yields the boundary condition 



•yr 



(A5) 



with the matching matrix 



V 



r de i(Tz<t>/2 
de~ iaA/2 f 



(A6) 



where = Xj?(0) — Xl(0). The matrix V satisfies the unitarity condition VV^ = 1, 
provided by the relations among the normal electron scattering amplitudes in Eq. (]A3): 



r = -d(r /<£)*, \d\ 2 = \d\ 2 = D, \r\ 



1^12 



R=l-D. 
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APPENDIX B: BOUND STATE CURRENT 



Equation ( |4.10| ) for the current of a single bound state can be derived directly |K| from 



the Bogoliubov-de Gennes equation Eq. ( [2.5|) , (j2.2|). The derivation is valid for junctions 
with an arbitrary non-superconducting region between superconducting electrodes. We as- 
sume for simplicity that the phase of the order parameter, Eq. ( |2.3| ), in the electrodes is 
constant and equal to ±0/2 in the right and left electrodes respectively. Let \l/(r, E) to be 
a normalized wave function of the Andreev bound state with energy E, 

H^-E^ = 0, ^(x = ±oo) = 0. (Bl) 

The energy and the wave function of the bound state depend on the phase difference <fi. 
Taking the derivative with respect to <fi of Eq. ([Bl]) and considering a scalar product of the 
resulting equation with the function we have 

Jd s r(*,±(H-E)*) = 0, (B2) 

where the brackets denote a scalar product in electron-hole space, similar to Eq. ( |3.8|) . In 
this equation the derivative of the Hamiltonian has the form 

dH dA zsigns ~ 

-dj = ^ = 2 G (B3) 

in accordance with Eqs. ( |2.2| ), (|2.3|). Substituting relation ( [B3| ) into Eq. ( [52]) and taking 
into account that the function \1/ is normalized, we get 

The continuity equation for the charge current, 

I(x, E) = ^-(p x - p' x ) J d 2 r ± ^(x'), *(x)) x=x ,, (B5) 



in accordance with Eq. ( |B1| ) has the form 

d 



i— I(x, E) = eJ d 2 r ± (*(x), [a z , A]V(x)). (B6) 
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Substituting relation ( p3|) into Eq. (|B6| ) and performing integration of this equation over 
the whole x axis we get 

1(0, E) = 2e J d 3 r(^,^^). (B7) 

In equation ( JB7| the current at infinity drops out because of decay of the bound state wave 
function, J(±oo) = 0. The magnitude 1(0) is formally taken in the middle of the junction; 
however, the current has the same magnitude along the whole non-superconducting region 
according to the conservation equation ( |B6|) . Comparison of Eqs. ( [Bp) and ([B7|) finally 
yields 

dE 

1(E) = 2e^. (B8) 
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FIGURES 



FIG. 1. SIS tunnel constriction. 

FIG. 2. Quasiparticle spectrum and position of the incoming states: 1(3) - hole (electron)-like 
quasiparticle incident from left; 2(4) - hole (electron)-like quasiparticle incident from right. 

FIG. 3. Spatial configuration of the edges of the superconducting energy bands in a long con- 
striction: E m i n , E max = ±A + p s (x)v. A potential well appears in upper (lower) band for electrons 
moving in a direction opposite to (along) the supercurrent. 

FIG. 4. Scattering diagram of voltage biased superconducting tunnel junctions. Solid (dotted) 
arrows indicate scattering in the normal (Andreev) channel. Filled triangles indicate supercon- 
ducting bound states. Transmission (reflection) occurs into side bands with odd (even) indices. 

FIG. 5. Scattering diagrams of voltage biased normal tunnel junctions, a) scattering of normal 
holes with spectrum E^ = fi — (p 2 /2m); represents an elementary fragment of the diagram in Fig. 
4 for j = 1. b) scattering of the normal electrons with spectrum E e = p 2 /2m. c) conventional 
diagram of elastic electron scattering in biased tunnel junctions; the local chemical potentials in 
the electrodes are then shifted by eV 

FIG. 6. Three kinds of processes contributing to the tunnel current at large bias eV S> A: a) 
creation of a real excitation across the gap by forward scattering; b) excitation of the Andreev bound 
state due to creation of a real excitation via backward scattering (dashed arrow); c) imbalance of 
ground state modes due to creation of a real excitation via backward scattering. Excess current is 
caused by processes b) and c). 
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FIG. 7. Scattering processes contributing to the subgap current, a) single particle scattering 
into side band n = 1 gives the main contribution at eV > 2A. b) excitation of the Andreev bound 
state (n = 1) due to backward scattering into side band n = 2 gives the main contribution at 
eV > A. c) single particle scattering into side band n = 3 and simultaneous excitation of the 
Andreev bound state in side band n = 1 gives the main contribution at eV > 2 A/3. 

FIG. 8. a) Density of states v{E) = \E n /£ >n \ of the side bands Eq,E\,E2 at applied voltage 
V > V\ (right), position of singularities of side band density of states as function of applied voltage 
for current I2 (left), ^ : E\ = ±A, 2~ : E2 = —A. b) Density of states of the side bands 
Eq,Ei,E2, £"3 at applied voltage V2 < V < V\ (right), position of singularities of side band density 
of states as function of applied voltage for current ^3 (left), 1^ : E\ = ±A, 2^ : E2 = ±A, 
3- : E 3 = -A. 



FIG. 9. Schematic picture of partial I n -V characteristics. 
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